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Abstract: In this paper, we present a new method for solving singularly perturbed two-point singular 
boundary value problems. Its exact solution is represented in a form of series in reproducing 
kernel space. In the mean time, the n-term approximation tn (x) to the exact solution u(x) 
is obtained and is proved to converge to the exact solution. Some numerical examples are 
studied to demonstrate the accuracy of the new method. Results indicate that the new 
method is simple and effective. 
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1 Introduction 


In this paper, we consider the following singularly perturbed two-point singular boundary 
value problem in reproducing kernel space 


H E u(x) = f(x x 
cul(a) + Powe) + ule) = f) 0<2S1, 


u(0) — 0, (1) 
u(1) = 0, 


where 0<e <1, u(x) € W2(0, 1], p(x) = O(z?) (a > 0) and q(x) = O(2°) (8 > 0). 

The singularly perturbed differential equations arise in a variety of differential applied math- 
ematics, fluid mechanics, quantum mechanics, optimal control, gas dynamics, nuclear physics, 
chemical reaction, studies of atomic structures and atomic calculations. Therefore, the prob- 
lem has attracted much attention and has been studied by many authors. In general, classical 
numerical methods fail to produce good approximations for these equations. Hence one has 
to go for non-classical method. Some non-classical methods have been suggested by various 
authors!!~4], But only few authors have developed numerical methods for singularly perturbed 
two-point singular boundary value problem. Mohanty and his co-workers have given the numeri- 
cal solution of singularly perturbed two-point singular boundary value problem using convergent 


tension spline method and non-uniform mesh tension spline method!®:41, 
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In this paper, we will give a representation of exact solution to (1) and approximate solution 
in the reproducing kernel space under the assumption that the solution to (1) is unique. The 
approach is simple and effective. 
For convenience, we take p(x) = z and q(z) = z? in (1). After multiplying (1) by x7, we 
find that 
ea?" (x) + zu'(z) + u(x) = z? f(z) O<e<, 
u(0) = 0, (2) 
u(1) — 0, 
where 0<e <1, u € W3[0, 1], z?f(z) € Wz[0, 1]. 
Clearly, the solution to (2) is the solution to (1). So we only need to obtain the solution of 
(2). Write F(z) = z? f(z) simply and put Lu = z?u" (x) + xu'(z) + u(x). Then (1) can further 
be converted into the following form: 
Lu(z) = F(z), OS zl, 
u(0) = 0, (3) 
u(1) = 0, 

where 0 «e <1, u € W3[0,1, F(z) € Wj[0, 1). W2[0, 1] and W3(0, 1] are defined in the 


following section. 


2 Several reproducing kernel spaces 


2.1 The reproducing kernel space W3[0, 1] 
The inner product space W3[0, 1] is defined as W3(0,1) = {u(z) | u, u', u" are absolutely 
continuous real value functions, u, u’, v", u®) € L?[0,1], u(0) = 0, u(1) = 0}. The inner 


product in W2[0, 1] is given by 


1 
(u(y),v(y)) ws = J (36uv + 49u'v' + 14u"v" + uy )dy, (4) 
3 0 


and the norm ||ul|ws is denoted by lullws = q4/(u,u)wg, where u, v € W3(0, 1]. 
Theorem 2.1 The space W3(0, 1] is a reproducing kernel space. That is, for any u(y) € 


W3D, 1] and each fixed x € [0, 1}, there exists Re(y) € W3(0,1], y € (0, 1], such that 
(u(y), Ra 9)); = uz). 
The reproducing kernel Rz(y) can be denoted by 


cie + cae + cae?! + cse?” + ege?! + cge*¥, — y X m, 
Rz(y) = (5) 
die + doe + dae?! + dye7?¥ + dse? + dee ?Y, y >=. 


The coefficients of the reproducing kernel R;(y) and the proof of Theorem 2.1 are given in [8]. 
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2.2 The reproducing kernel space W}(0, 1] 

The inner product space W1[0, 1] is defined by W1[0, 1] = (u(z) | u is absolutely continuous 
real value function, u, u’ € L?[0,1]}. The inner product and norm in W2[0,1] are given, 
respectively, by 


1 
(vr), = | (uim, dla = eng 


where u(x), v(x) € W3[0,1]. In [7], the authors have proved that W2[0,1] is a complete 
reproducing kernel space and its reproducing kernel is 


R.(y) = saan [cosh(z + y — 1) + cosh (|x — y| — 1) ]. 


3 The solution of (3) 


In this section, the solution of (3) is given in the reproducing kernel space W310, 1]. 

In (3), it is clear that L : W3[0,1] —^ W2{0, 1] is a bounded linear operator. Put y;(z) = 
Rz,(x) and vi(z) = L*p;(x) where L* is the adjoint operator of L. The orthonormal system 
[v;(z) 2, of W3[0, 1] can be derived from the Gram-Schmidt orthogonalization process of 
{bi(x) 21 


i 
Vir) = Y Binve(z), Bix >0, i—L2,-. (6) 
k=1 
Theorem 3.1 For (3), if {z;}9, is dense on [0, 1], then (;(x));2, is the complete system 
of W3[0, 1] and (x) = L,Rs(y)ly-z;- 
Proof We have 


pila) = (Lyi) (2) = (Q9) (v), R«(y)) = (vily), Ly Re(y)) = LyRe(y)ly=as- 


The subscript y by the operator L indicates that the operator L applies to the function of y. 
Clearly, ;(z) € W3[0, 1]. For each fixed u(x) € W32[0, 1], let (u(z), v;(z)) 20, i = 1,2,---, 
which means that 


(u(z), (^5) (2)) = (Lu) e) Qa) (2) = 0. (2) 


Note that {2;}S2, is dense on [0,1]. Hence, (Lu)(x) = 0. It follows that u = 0 from the 
existence of L71. So the proof of Theorem 3.1 is completed. 

Theorem 3.2 If (z;)j2, is dense on [0,1] and the solution of (3) is unique, then the 
solution of (3) is 


u(x) = bD 3 Bir f (xe) bi(x). (8) 
i=1 k=1 


Proof Applying Theorem 3.1, it is easy to know that (V;(z))2, is the complete orthonor- 
mal basis of W3(0, 1]. 
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Noting that (v(x), yi(x)) = v(z;) for each v(x) € W3[0, 1], we have 





ulz) = Y ^ (u(z), Yi) Bie) = Y. Y. Ba (ule), Lex (2) 3, (x) 
i=l i=1 k=1 
= SO Y Bik (Luz), ex(2)) 9; (x) = Y, Y. Ba (F), prle): n) 
t=1 k=1 i=1 k=1 
= M 6af(x)9 (x). (9) 
i=1 k=1 


So the proof of Theorem 3.2 is completed. 
Now, the approximate solution u,(x) can be obtained by the n-term intercept of the exact 


solution u(x) and 


Un (z) = Y Y Ba f (zk)v; (x). (10) 


i=1 k=1 
Theorem 3.3 Assume u(x) is the solution of (3) and r„(£) is the error between the 
approximate u,(x) and the exact solution u(x). Then the error rp (x) is monotone decreasing 


in the sense of || - || ws. 
Proof From (9), (10), it follows that 


ls lg = | 95 XS Anf) lu; 
i=nt+1 k=1 
oo i 2 
= M (bese). (11) 
i=n+1 k=1 
(11) shows that the error r, is monotone decreasing in the sense of ll - llwg- So the proof of 


Theorem 3.3 is completed. 


4 Numerical example 


In this section, some numerical examples are studied to demonstrate the accuracy of the 
presented method. The examples are computed using Mathematica 4.2. Results obtained by 
the method are compared with the exact solution of each example and are found to be in good 
agreement with each other. 

Example 1 Consider the iyd 


eu” (x) + zu "(z) + au) = f(z), O<2<1, 


u(0) = me 0, 
u(1) — 0, 
where f(x) = 2 — 2e — 3. The true solution is x — z?. Using our method, we choose 26 points 


on [0, 1] and take € = 107°, 107, 10779, respectively. The numerical results are given in Table 
1, Table 2 and Table 3. 
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Table 1: Numerical results for Example 1 (n = 26, € = 107?) 










absolute error 










0.2304 0.2306 
0.1600 0.1600 
0.0384 0.0383996 





Table 2: Numerical results for Example 1 (n = 26, € = 1079) 








absolute error 


| sic 





Table 3: Numerical results for Example 1 (n = 26, € = 1071?) 






absolute error 
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Example 2 Consider the equation 








1 
" / MET = <1 
ew" (x) + cm (z) + zz ule) f(x), O<a2<1, 
u(0) = 0, 
u(1) = 0, 
where 
oe reos) CSCE 02 dub 2 


The true solution is sin(zz). Using our method, we choose 26 points on [0,1] and take € = 
1072, 1076, 10-15, respectively. The numerical results are given in the following Table 4, Table 
5 and Table 6. 


Table 4: Numerical results for Example 2 (n = 26, € = 10^?) 











T approximate solution u25 absolute error 
0.0001 0.0314157 1.8E-09 
0.08 0.31868 46805 
0.16 0.481754 0.481702 5.1E-05 
0.32 0.844328 0.844275 5.2E-05 





0.48 0.998027 0.997974 5.2E-05 
0.64 0.904827 0.904776 5.1E-05 


0.80 0.587785 0.587736 4.9E-05 
0.96 0.125333 0.125274 5.9E-05 
0 


1.00 0 0 











Table 5: Numerical results for Example 2 (n = 26, € = 10~°) 















absolute error 





true solution u(x) approximate solution u25 


0.998027 0.997974 


0.125333 0.125274 
0 0 
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Table 6: Numerical results for Example 2 (n = 26, e = 10 19) 
true solution u(x) approximate solution ues absolute error 
0.0001 0.0314159 0.0314157 1.8E-09 
d een a 
0.08 0.24869 0.248643 4.6E-05 
0.481754 0.481702 5.1E-05 
0.32 0.844328 0.844275 5.2E-05 
0.998027 0.997974 5.2E-05 
0.904827 0.904776 5.1E-05 
0.587785 0.587736 4.9E-05 
0.125333 0.125274 5.9E-05 
1.00 0 0 0 
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MEFRES 1 75 A 
*x deo 


(1- PARE TK FE BCS IR, M/K 150001; 2- MAME AERP, PR 210045) 


iff 


EERIE, RIMA- ARRA REED PS OZ RE Z3 S. "CURRERE TE EAE E i nf E 


EAZCBOE XUI. INS, FE nME us (x) JE T8 8858. u(x), FF AAT ALE us (x) KAF u(x). R 
ATE UES HH EB US UE TEE, FR A, Sha T ARA AA 
Kid: HA: ARROW: BER 


